arXiv: 1502.03605v2 [astro-ph.EP] 24 Apr 2015 


Astronomy & Astrophysics manuscript no. paper_domETAL_revised2 ©ESO 2015 

April 27, 2015 


Can we constrain interior structure of rocky exoplanets from 

mass and radius measurements? 

Caroline Dorn 1 , Amir Khan 2 , Kevin Heng 1 , Yann Alibert 1 , James A. D. Connolly 3 , Willy Benz 1 , and Paul Tackley 2 

1 Physics Institute, University of Bern, Sidlerstrasse 5, CH-3012, Bern, Switzerland 
e-mail: Caroline . dorngspace . unibe . ch 

2 Institute of Geophysics, ETH Zurich, Sonneggstrasse 5, 8092 Zurich 

3 Institute of Geochemistry and Petrology, ETH Zurich, Clausiusstrasse 25, 8092 Zurich 


Received September 2014 


ABSTRACT 

Aims. We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the 
form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial 
exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. 

Methods. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain 
confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated 
uncertainties. 

Results. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent 
estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to 
exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision 
required to characterize the interior of exoplanets. 

Conclusions. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental 
abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) 
the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement 
accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abun¬ 
dances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. We provide a general 
methodology of analyzing interior structures of exoplanets that may help to understand how interior models are distributed among 
star systems. The methodology we propose is sufficiently general to allow its future extension to more complex internal structures 
including hydrogen- and water-rich exoplanets. 
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1. Introduction 

Major advances in detection and characterization of exoplan¬ 
ets have been achieved over the past decade. To date, several 
hundred have been characterized in terms of mass and radius by 
space-based or ground-based observations. Continued improve¬ 
ments in observational techniques allow for more precise infer¬ 
ence of mass and radius. To date, even a few small-mass exo¬ 
planets (< 10 Earth masses (Me)) with uncertainties in mass 
and radius below 20% and 10%, respectively, have been dis¬ 
covered (e.g., Kepler-36b, CoRot-7b) (e.g., Pepe et aT7| 201 3| >. 
Space-based missions, like Kepler (Borucki et al. 2010| > and 
CoRoT (Borde et al. 2003), have been able to measure the tran¬ 
sits of these objects and infer the radius. Within the first 16 
months of the Kepler mission alone, 207 planetary candidates 
of radius R < 1.25 Earth radii (Re) and 680 super-Earth-sized 
(1.25i?s < R < 2 Re) planetary candidates were reported by 
|Batalha et al.| ( |20T3] >. 

From knowledge of mass and radius we are able to derive 
constraints on the interior structure of exoplanets. Previous stud¬ 
ies concerned with rocky exoplanets (e.g., Valencia et al.[2006; 


computing mass-radius relations based on terrestrial-type inte¬ 
rior structures and compositions. These studies showed that dif¬ 
ferent interior models are capable of explaining the observations 
within their uncertainties. Such “forward" approaches, however, 
do not quantify the inherent degeneracy of interior structure 
models. For example, Valencia et al.| (j2006 2007c)i; |Sotin et ah] 
(2007); Wagner et al. ( |2011[ > it has been shown that different 
core sizes and mantle compositions affect the mass-radius rela¬ 
tionship, but to our knowledge no comprehensive study of the 
degeneracy of these model parameters has been performed. As a 
consequence, it is not fully understood which interior-structure 
parameters can be estimated from the data and which structure 
parameters tend to be strongly correlated. 

In the light of the inherent ambiguity in determining interior 
structure from mass and radius (e.g., |Howe et al.|2014||Rogers &| 
Seager|2010|i, analysis is most sensibly conducted in an inverse 


Sotin et al.|2007[|Seager et al.|2007[|Fortney et al.|2007H Wagner 

, Mocc 


et al. 2011 1 Mocquet et al,|2014|> have generally concentrated on 


sense ( [Mosegaard et al.|2002| . We therefore propose a complete 
Bayesian inverse analysis by employing a Markov chain Monte 
Carlo (McMC) method to provide full probability distributions 
for the model parameters of interest (mantle composition and 
core size). 
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Early studies of mass-radius relations ( |Zapolsky & Salpeter| 
1969) considered only homogeneous and simple monoatomic 


bulk compositions (H, He, C, Fe, Mg, and H/He mixtures). More 
recent studies assume differentiated terrestrial-like compositions 


(e.g., Valencia et al. 

2006 Sotin et al. 2007; Seager et al. 2007 ; 

Fortney et al. 2007 

. Based on solar system observations and 


planet formation and differentiation models, terrestrial planets 
are generally thought to be differentiated in an iron-rich core, a 
silicate mantle and a crust (see [Howe et al.|2014[ for a review of 
previous and current work). In this work, we focus on terrestrial- 
type rocky planets and assume that these consist of a pure iron 
core, a silicate mantle comprising the oxides Na 20 -Ca 0 -Fe 0 - 
Mg 0 -Al 2 03 -Si 02 , and that volatiles have a negligible effect on 
mass and radius. 

In order to be able to compute the radius for a planet of a 
given mass and composition (the ’’forward" model) we solve 
the hydrostatic equilibrium equation coupled with a thermody¬ 
namic approach based on Gibbs free-energy minimization and 
Equation-of-State (EoS) modeling. In this regard, our work is 
similar to earlier studies ( Valencia et al.|2007b||Sotin et al.|20(T7| 


Seager et al. 2007[|Fortney et al.|2007 1 , except for the use of ther¬ 
modynamic modeling. This addition allows us to compute from 
first principles of thermodynamics, the density profile of the 
planet for an arbitrary bulk silicate mantle composition. More¬ 
over, given that cosmochemically-derived elemental abundances 
(e.g., Fe, Mg, and Si) have proved important as constraints on the 
composition of the terrestrial planets (e.g., Grasset et al.] |2009[ 
Rogers & Seager] [2010[ >, we investigate the effect of including 
abundance constraints in determining interior structure. 

Although similarities between the present work and the 
Bayesian analysis of| Rogers & Seager[ ( |2010j > exist, our approach 
nonetheless differs in several respects: 


- our inversion method is applicable to planets of any mantle 
composition and core structure; 

- our inversion method is general in that it can be easily 
adapted to more complex planetary models that include 
oceans and atmospheres; 

- we infer probability distributions on interior-structure pa¬ 

rameters (core radius and mantle composition) for a wide 
range of masses (O.IM^ < M < 10 Me) and radii 

(0.5Re < R < 1.7 Re)', 

- we demonstrate how observations of a few key elemen¬ 
tal abundances (e.g., Fe, Mg, Si) of the host star’s photo¬ 
sphere significantly reduces the degeneracy of interior struc¬ 
ture models; 

- we quantify the influence of measurement uncertainties on 
the determination of interior structure models. 


The outline of this study is as follows: In Section[2]we intro¬ 
duce the physical model that links data and model parameters, 
and outline the inversion strategy. In Section [3] we first test our 
method on Earth using Earth bulk elemental abundances derived 
from geochemical analysis of mantle rocks. Next, we apply the 
method to all terrestrial solar system planets. In Section [4] we 
present results for synthetic exoplanets and Kepler-36b followed 
by discussion and conclusions. 


2. Methodology 

2.1. Abundance constraints 

In summary, the importance of additional constraints in the 
form of chemical abundances of the elements Fe, Mg, Si are 
investigated in detail. Arguments from planet formation favor 
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certain interior compositions as discussed by e.g., |Grasset et al.| 
( |2Q09| >; Rogers & Seager ( 2010| ). 

Observations and theoretical considerations suggest that exo¬ 
planet bulk composition is dictated by stellar composition: First, 
relative elemental abundances of Fe, Mg, and Si are similar 
among the Sun, Earth, Mars, the Moon as well as meteorites 
(|Lodders||2003| |Drake & Righter |2002| |McDonough & S~un| 
1995|l. Meteorites are believed to be chemically similar to the 


building blocks of planets (Morgan & Anders 1980] >, both be¬ 
ing condensates from the solar nebula that experienced the same 
fractionation processes. Second, it has been demonstrated that 
during planet formation, planetary bulk and stellar ratios of Fe/Si 
and Mg/Si are very similar (e.g., |Bond et al.| 2010[ |Elser et ST] 


|2012[ [Johnson et al.|2012[|Thiabaud et al.|2014[ > since the refrac¬ 

tory elements Fe, Si, and Mg condense at similar temperatures 
corresponding to small distances to the host star (~ 1 AU). 

For the Sun and our solar system planets, bulk Fe/Si and 
Mg/Si (henceforth Fe/SiQ and Mg/SiQ) have been inferred 
from photospheric analysis and from geochemical analysis of 
meteoritic and cosmochemical material, respectively. For extra¬ 
solar systems, we shall assume that Fe/SiQ and Mg/SiQ equal 
to their stellar counterparts Fe/SL*. and Mg/SL*.. 


2.2. Data 

The aim of our study is to characterize the interior structure 
of terrestrial exoplanets based on inferred mass and radius from 
radial velocity and transit observations and constraints on bulk 
Fe, Mg, and Si abundances (Fe/SiQ and Mg/SiQ) inferred from 
stellar photospheric analysis. Mass, radius and bulk abundance 
constraints constitute our data d. For the purpose of illustrating 
the relative importance of data on inverted model parameters, we 
consider the following two cases: 

Case A Data are mass M and radius R. 

Case B Data are mass M, radius R , and compositional con¬ 
straints (Fe/SiQ and Mg/SiQ) obtained from stellar spec¬ 
troscopic observations of the host star. 

Note that Fe/SiQ expresses the mass ratio between the mass 
of iron to silicate for the entire planet (core and mantle). Given 
that we assumed a pure iron core, this ratio is equal to the mass 
of iron in core and mantle divided by the mass of silicate in the 
mantle M s i,mantle- This can be written as follows: 


Fe/SiQ = Fe/SiQ + M core /M s - h mantle (1) 

where M core is mass of iron in the core, and Fe/SiQ denotes 

the mass ratio of Fe/Si in the mantle. Mantle Mg/Si is referred to 
as Mg/SiQ. Since all Si and Mg is assumed to be in the mantle, 

Mg/Si 0 = Mg/SiQ. 


2.3. Model parameterization 

Our exoplanet interior model consists of a layered sphere 
with an iron core surrounded by a silicate mantle (see Figure [TJ. 
Model parameters include core size R c and mantle silicate com¬ 
position using the NCFMAS model chemical system that com¬ 
prises the oxides Na 2 0 -Ca 0 -Fe 0 -Mg 0 -Al 2 03 -Si 0 2 and ac¬ 
counts for more than 98% of the mass of Earth’s mantle ( |Irifune| 
|1994] > (see Appendix |C|. We parameterize mantle composition 
by mantle Si-content defined as Msi.mantie/A^mantie. (referred 
to as SIq), Fe/SiQ and Mg/SiQ. Na 2 0, CaO, and AI 2 O 3 
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^surf, ^surf 



Fig. 1. Illustration of model parameterization. Parameters are core radius R c and mantle composition c comprising the oxides Na20-Ca0-Fe0- 
Mg0-Al 2 03-Si02- Mantle temperature profile T, surface temperature T sur f and pressure P sur f are fixed input parameters. In the iron core an 
adiabatic temperature profile is used. 



Fig. 2. Inversion scheme. Note that the bulk compositional constraints from the host star (case B) used both in prior and data are independent. The 
iterative scheme is repeated until the output converges to the posterior distribution. See section[2]for more details. 


are minor components and their relative abundance can be con¬ 
strained from stellar photospheric observations (e.g., |Elser et al.| 
2012). In this study, we assume chondritic abundances (see Table 
TJ. The set of parameters that constitutes our model parameter m 
is thus: 

- core radius R c , 

- mantle Si-content Sigj, 

- mantle Fe/SiQ, 

- mantle Mg/SiQ. 

The mantle composition parameters Si^, Fe/Sigj, and 
Mg/SiQ depend on each other, i.e. the fractions of oxides in 
the NCFMAS system must sum up to one. All parameters are 


assumed to be uniformly distributed unless stated otherwise (see 
Table [Tj). 

2.4. Prior information 


For case A, prior bounds have been chosen based on: (1) 
the assumption that the bulk metallic composition of a planet is 
similar to the abundance of the host star ( jThiabaud et al.|2014| l 
and (2) the fact that chemical abundances of observed stars with 
planetary companions seem to fall within certain ranges (e.g., 
|Beirao et al.|2005[|Gilli et al.|2006| l. Thus, the variability of ob¬ 
served elemental abundances in stars defines the prior bounds on 
Fe/Si@ and Mg/Si@. 
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Gilli et al. 2006 


Table 1 . Prior ranges of model parameters for case A (Beirao et al. 2005 
Uniform distributions unless stated otherwise. Si, Fe, and Mg are in wt% (weight percentage). Minor components of Na 20 . 
fixed to chondritic abundances of 2 wt%, 2 wt% and 4.8 wt%, respectively [Lodders|2003[ . 


Thiaba ud et al.|20~l4 1 and case B |Lodders|2003| >. 

CaO, and AI 2 O 3 are 


model parameter of m 

case A 

case B 

Si ® 

8 -34 (non-uniform^) 

14 -30 (non-uniform') 

Fe/SiQ 

0-4.15 

0 -2.23 

M g/ Si (o) 

0.61-2.32 

0.89 ± 0.08 (Gaussian) 

core radius R c 

0 -R 

0 -R 

f see appendix [a] 


For case B, prior bounds on Fe/SiQ and Mg/SiQ are de¬ 
rived from Fe/SiQ and Mg/SiQ inferred from the host star’s 
photospheric abundance: (1) Since all Si and Mg is assumed to 
be in the mantle, Mg/SiQ defines the prior bounds on Mg/SiQ; 

(2) Fe, on the other hand, is distributed between core and man¬ 
tle. Thus, the bulk constraint Fe/SiQ defines only the upper 
bound of the prior on Fe/SiQ. For exoplanets, Fe/SiQ and 

Mg/SiQ are defined by the host star’s photospheric elemen¬ 
tal abundance since we assume that Fe/SiQ = Fe/SL*. and 
Mg/SiQ = Mg/SL*.. Mg/SL*. is Gaussian-distributed. We con¬ 
sider solar estimates for Fe, Si, and Mg ( |Lodders12003) ) and as¬ 
sociated uncertainties as representative of Fe/SiQ and Mg/SiQ 
throughout this study unless stated otherwise. 

In both cases, the prior bounds on SiQ depend on the prior 
bounds of Fe/SiQ and Mg/SiQ, since the fractions of oxides 
in the NCFMAS chemical system must sum to one. This implies 
a non-uniform prior of S1 q (see Appendix |A|>. 

Additional input parameters that are fixed in this study are 
( 1 ) surface pressure and temperature, ( 2 ) thermal gradient in the 
mantle based on adiabatic gradient of Earth’s mantle, and (3) 
absolute amount of minor components (Na 20 , CaO, and AI 2 O 3 ) 
in the mantle. 


Using a fixed temperature profile for the mantle is demon¬ 
strably not a limitation of our model. Many authors have 
shown that the variations due to temperature for rocky planets 
are small compared to uncertainties of material properties of 
silicate or iron compositions (e. g., |Sotin et al.| 2007| Seager| 
|et al.||2007] |Grasset et al.||2009| [Valencia et al.||2010| i' This is 
demonstrated in Appendix [B] where it is shown that varying 
mantle temperatures introduces little variation in computed 
densities. As a consequence, we fix temperature profiles by 
using Earth’s surface temperature, lithosphere and mantle 
temperature gradient. 


Equation-of-State (EoS) modeling. The complete solution of the 
forward problem is summarized as follows: 

{c,R cl T} a -\M?\ p^{M,R} (3) 

where c denotes the NCFMAS mantle composition, T tempera¬ 
ture, A4 equilibrium mantle mineralogy, and p density (all pa¬ 
rameters are a functions of radius). The forward operator gi 
embodies Gibbs free-energy minimization that computes stable 
mineral modes as a function of P (pressure), T, and composi¬ 
tion, g 2 calculates density, g 3 integrates p over R to compute M. 
We assume subsolidus conditions throughout. 


2.5.1. Thermodynamic modeling 


Possible mantle compositions are explored within the Na 20 - 
Ca 0 -Fe 0 -Mg 0 -Al 203 -Si 02 model system. To compute the 
mantle density profile, we employ a self-consistent thermody¬ 
namic method. For this purpose we assume: thermodynamic 
equilibrium; mantle mineral phases of large exoplanets are those 
that potentially occur in the deep mantle of the Earth; and, when 
required, the thermodynamic properties of the mineralogy can 
be extrapolated to more extreme conditions than realized on 
Earth. The equilibrium mineralogy and, consequently its den¬ 
sity, is computed as a function of pressure, temperature, and bulk 
composition by Gibbs energy minimization ( jConnolly||2009| l. 
For these calculations the pressure is obtained by integrating 
the load from the surface boundary and temperature is obtained 
by integrating an Earth-like temperature gradient from an arbi¬ 
trarily assumed surface temperature. Mineral equations of state 
and parameters are as given by Stixrude & Lithgow-Bertelloni 
( |2005l[20TT) . Because the equilibrium assumption is dubious at 
low temperature (e.g.. Wood & Holloway|1984| , for models that 
require rock properties at temperatures below 800 K, the stable 
mineralogy is first calculated at 800 K and its physical properties 
are then computed at the temperature of interest. 


2.5. Interior structure model 


2.5.2. EoS for iron core 


Data and model parameters are linked by a physical model 
embodied by the forward operator g( - ) 


d = g(m) 


( 2 ) 


For a given model m, interior structure (density profile), the 
total mass, and Fe/SiQ are computed for the purpose of com¬ 
paring with observed data d. Equation [^represents the forward 
problem, which is computed using thermodynamic method and 


In order to compute the core density profile, we use the EoS 
for pure iron derived by Belonoshko ( 2010|), which i s similar to 
the Mie-Griineisen-Debye EoS ( | Jack son & Rig den 1996), ex¬ 
cept for a different thermal pressure term. This EoS for iron 
is preferred, as it allows us to closely reproduce high-precision 
volumetric experiments of iron under high-pressure and high- 
temperature conditions. 

The adiabatic gradient in the core is computed as 


dT _ jT 


(4) 


Article number, page 4 of 


19 












































Dorn et al.: Can we constrain interior structure of rocky exoplanets from mass and radius? 


where Kg is adiabatic bulk modulus, which is related to the 
isothermal bulk modulus Kt through Kg = (1 + 7 aT) Kt- 
The parameters of the Mie-Griineisen-Debye equation 7 = 

7o f° r i ron are Po = 8.334 g/cm 3 , q = 0.489, and 

70 = 2.434. Density p and thermal expansion coefficient a are 
then directly computed from Belonoshko’s EoS. 


2.6. Inversion method 


Probabilistic inversion methods are suitable when data 
are sparse, the physical model (Eq. [2]i is highly non-linear 
and/or it is expected that very different models are consis¬ 
tent with information and data (e.g., Mosegaard & Tarantola 
1995 ). While stochastic sampling-based approaches for high¬ 
dimensional problems are computationally expensive, numerous 
global search methods exist. These include Markov chain Monte 
Carlo methods (McMC), simulated annealing or genetic algo¬ 
rithms. McMC has the advantage that parameter uncertainties 
are formally assessed. Furthermore, the McMC method gener¬ 
ates samples whose distribution converges to a stationary dis¬ 
tribution coinciding with the posterior probability density func¬ 
tion (pdf) (e.g., Hastings |l970j ). Although the dimensionality of 
the presented inverse problem is low (4 parameters), offering 
the possibility of performing grid-search, we nonetheless em¬ 
ploy McMC, because we intend to expand the complexity of the 
problem to higher dimensions, in the future. For consistency, we 
have also used the grid-search method to verify our results. 


2.6.1. Bayesian analysis 

We employ a Bayesian method to compute the posterior 
probability density function (pdf) for each model parameter (R c , 
Si 0 , Fe/Sigj, Mg/Sigp from data (M, R and Fe/Sio) and 
prior information. According to Bayes’ theorem, the posterior 
distribution for a fixed model parameterization m, conditional 
on data d, is given by (Tarantola & Valett e|T982 ): 


p(m|d) oc p(m)L(m|d), (5) 

where p(m) represents prior information on model parame¬ 
ter m. L(m|d) is the likelihood function and can be interpreted 
in probabilistic terms as a measure of how well a model m fits 
data d. 

Assuming uncorrelated and normally-distributed residuals 
(see Eq. [7j, the likelihood defines a “measure” between d and 
m and can be written as: 


L(m|d) = 


(2n) N /2(nh*?) 1/2 


exp - 



( 6 ) 


2.6.2. Markov chain Monte Carlo (McMC) 

For high-dimensional and non-linear inverse problems or 
sparse data, it is in practice impossible to derive the posterior 
distribution analytically. McMC techniques offer an efficient 
method of performing Bayesian analysis (Eq. [5]). McMC sam¬ 
pling methods iteratively search the space of feasible solutions. 
The iteration scheme used is based on the Metropolis algorithm, 
which proceeds as follows: 

1. An initial starting model m o]( | is drawn at random by sam¬ 
pling from the prior distribution. 

2. The posterior density of m 0 id is calculated by evaluating the 
product of the likelihood of the corresponding forward model 
and prior density. 

3. The current model parameter is perturbed to obtain m new , 
which is subsequently created from a proposal distribution. 
Here, the proposal distribution is uniformly bounded and 
centered around mold- Generally, it is chosen such that the 
size of the model perturbation allows for a reasonable rate 
of accepted transitions in the McMC procedure, typically 
around 30% ( |Gilks et al.|1996| l. 

4. The proposal m new is accepted with probability (Mosegaard 
|& Tarantola|1995] l: 

Paccept — rnin{l, exp(Z(m new |d) /(m 0 id|cl))j' 

5. If the proposal is accepted the Markov chain moves to m new , 
otherwise the chain remains at m 0 id. 

6 . Steps 2 to 5 are repeated. 

After many iterations, the accepted samples that are gener¬ 
ated with this approach are ( 1 ) independent of the starting model 
and (2) distributed according to the posterior distribution. The 
efficiency of sampling is strongly dependent upon the proposal 
distribution. If this distribution is badly chosen, the acceptance 
rate of solutions might be unacceptably low, resulting in very 
poor efficiency regarding convergence (Hastings 1970[ l. On the 
contrary, if the proposal distribution is well chosen, the McMC 
sampler will rapidly explore the posterior target distribution. 
Here, the posterior information is gathered from a large num¬ 
ber of sampled models (~ 10 4 ). Only the statistical nature of 
sampled model features are of interest. Presently, we will con¬ 
centrate on computing marginal posterior distributions in order 
to provide the reader with a notion of model parameter uncer¬ 
tainties. Information on single parameters are obtained by one¬ 
dimensional (ID) marginals; higher dimensional marginal pdfs 
reveal the correlation that exists among several parameters. If 
the data are able to constrain the model (i.e. informative), differ¬ 
ences between prior and posterior pdfs are expected. 

Given the large number of models that have to be computed, 
the calculation of the forward model must be computationally 
very efficient. In our work, generating the planet’s internal struc¬ 
ture takes on average 0.8 seconds of CPU time on a four quad- 
core AMD Opteron 8380 CPU node with 32 GB of RAM. In all, 
we sampled about 10 ' models and retained around 10 4 models 
for further analysis. 


where N is total number of data points, rr, is the estimated error 
of the i th datum, and Q, is data misfit, or squared residual, given 
by 


Qi = (gi{ m) - di) T (g,;(m) - di) 


where g{ m) is the forward model discussed in section 2.5 


(7) 


3. Method validation 

3.1. Synthetic case 

We apply our method to a synthetic planet that ultimately 
serves as a ’’benchmark" test. Structure and composition of the 
synthetic model are defined as: 

- SiQ = 23.4 wt% 
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- Fe/SiQ = 0.5; Fe/Si 0 = 1-69 

- Mg/Si @ = Mg/Si 0 = 0.89 

- R c = 0.47 * R 


Using the interior structure model (Section [23| we calculate the 
synthetic data to be M = 1 Me and R = 1 Re- For the in¬ 
version we consider the following data uncertainties: ap = 0 , 
(7m = 0.001 M, together with Fe/SiQ = 1.69 ± 0.18 and 
Fe/SiQ = 0.89 ± 0.08. Results are shown in Figure [3] and 
demonstrate the ability of the method to retrieve the internal 
structure of the synthetic planet. As expected, the most likely 
posterior model overlaps with the true model. The inherent de¬ 
generacy in interior structure is apparent from the spread of the 
posterior. 


3.2. Solar System planets 

We test our method on the terrestrial solar system planets 
Earth, Moon, Mars, Venus, and Mercury, for which a series of in¬ 
dependent model parameter estimates are available from a host 
of geophysical, geo-, and cosmochemical data gathered in situ 
as well as from orbiting spacecraft (see Table [2]). These tests are 
important inasmuch as they allow us to benchmark our method. 
For simplicity, Table[2]only summarizes a single best estimate of 
independent interior models. Associated uncertainties are very 
different depending on data and model assumptions. For exam¬ 
ple, core radius R c of the Moon, Mars and Mercury, is estimated 
to be 310-320 km ( Kuskov & Kronrod|2001|i, 1680± 150 km 


( Khan & Connolly|2068| l, and 2020 ± 100 (Padovan et al. 2014| >, 

respectively. 

In the following, we extensively test different scenarios for 
the Earth: case A and B. Input for cases A and B is defined in 
section [2] These cases show that the compositional constraints 
Fe/SiQ and Mg/SiQ are crucial for determining interior struc¬ 
ture and that Fe/SiQ and Mg/SiQ are good proxies for them. 
How well data for case B allow us to recover the interior struc¬ 
ture is demonstrated for the Moon, Venus, Mars, and Mercury. 
For all terrestrial planets we use 1% and 0% uncertainties on 
observed M and R, respectively, and Fe/SiQ = 1.69 ± 0.18 
( |Lodders|2 003 ). Minor components of Na 20 , CaO, and AI 2 O 3 
are fixed to chondritic values (Lodders 2003). Prior ranges for 
the various parameters are listed in Table [T] 


Earth Results for the Earth for cases A and B are shown 
in Figure [4] Independent estimates of mantle composition and 
core size (Table [2]» are shown as green bars. Inversion of mass 
and radius only (case A) provides little information on mantle 
composition, but some constraints on R c . In contrast, addition 
of Fe/SiQ and Mg/SiQ (case B) enables us, as expected, to 
constrain mantle composition and thereby obtain improved es¬ 
timates of core size. Moreover, better agreement with the inde¬ 
pendent estimates is also observed. 

There is an underestimation of the core radius due to our 


assumed core composition. In the PREM model (Dziewonski 
|& Anderson| [l981 ), the core is less dense which reflects the 
presence of lighter elements (e.g., S, Si, C, O) in addition to 
Fe-Ni. Since we consider pure iron, we systematically overes¬ 
timate core density and therefore underestimate core size. For 
both cases A and B, core sizes are 0 < R c < 2930 km and 
1470 < R c < 3300 km (95% credible intervals), respectively. 
Furthermore, we observe that Mg/SiQ, is not well constrained 
by data. 
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Figure[5]shows the correlation between individual model pa¬ 
rameters for case A and B. In contrast to case A, strong corre¬ 
lations in case B are evident between model parameters: SiQ, 

Fe/SiQ, and R c . This is due to the constraint provided by 

Fe/SiQ that couples core and mantle, i.e., a larger core must 
be compensated by a lower Fe/SiQ or a larger SiQ. No signif¬ 
icant correlation of Mg/SiQ to the other model parameters is 
observed, because data are not able to constrain Mg/SiQ well. 


Moon Results for the Moon are shown in Figure [ 6 ] As for 
Earth, we see that SiQ, Fe/SiQ, and R c are well-constrained, 
i.e. prior and posterior pdf significantly differ, except for 
Mg/SiQ,. Predicted core size matches the independent estimate 
(Table [2ji reasonably well with 0 < R c < 380 km (95% credible 
interval), which incidentally also assumes a pure iron core. The 
mantle composition appears not to be fully retrieved here. This is 
because the Moon has the lowest Fe/SiQ among the terrestrial 
planets, including chondrites and satellites of the outer solar sys¬ 
tem ( [Kuskov & Kronrod||2001 1 . Therefore, the assumption that 
Fe/SiQ = Fe/SiQ leads us to overestimate Fe/SiQ and under¬ 
estimate SiQ. The independent estimate of mantle compositions 
nonetheless lies within the range of inferred posterior models. 


Mars Results for Mars are shown in Figure [7] Apart from 
Mg/SiQ, mantle composition (SiQ and Fe/SiQ) and core ra¬ 
dius (R c ) appear to be well-constrained. Although independent 
geophysical estimates of Khan & Connolly ( |2008[ ) lie within 
the posterior distribution, there is discrepancy that relates to 
the assumption of sulfur in the core model of Khan & Con- 
nolly ( 2008| >. As a consequence, we underestimate core size to 
0 < f? c < 1330 km (95% credible interval), which leads to an 
overestimate of Fe/SiQ and an underestimate of SiQ . Com¬ 


paring our results to models with pure iron cores (Rivoldini et al. 
|201 l[|Khan & Connolly|2008]>, reduces discrepancy. 


Detailed independent analyses of the mantle composi- 

Aitta: (|2012[>, e.g. impose 


Venus 

tion of Venus are not yet available, 
an Earth-like mantle composition and find that the composition 
of Venus’ core is likely enriched in lighter elements above that 
of Earth’s core to fit mass and radius. Consequently, as shown 
in Figure [ 8 ] we underestimate core radius to 0 < R c < 2810 
km (95% credible interval) compared to the independent esti¬ 
mates of Aitta|( |20T2) >. Generally, R c is well constrained, while 
the mantle composition (SiQ, Fe/SiQ and Mg/SiQ) is only 


moderately-to-weakly constrained. 


Mercury Mercury has a high bulk density and is thought to 
have experienced a different accretion scenario or undergone 
a different post-formation process in comparison to the other 
terrestrial planets ( |Cameron et al. 1988| |Strom & Spr ague 
2003). An early giant impact, is likely to have stripped off a 


large part of its silicate crust and mantle (Benz et ak] |1988| >, 
leaving behind a planet with a large iron core. Hence, the 
constraint Fe/SiQ = Fe/SiQ is not applicable to bulk Mercury 
as a consequence of which it is not possible to match data. 
We nonetheless show the posterior model parameter pdfs in 
Figure [9] but emphasize that these models have near-zero 
likelihoods. To improve fit to M a larger iron core is required. 
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Fig. 3. Sampled ID marginal posterior distributions for the synthetic test where model parameters are inverted for an Earth-sized planet: prior (red), 
posterior (blue) and independent estimates (green) of model parameters: (a) mantle Si-content SIq, (b) mantle Fe/SiQ, (c) mantle Mg/SiQ, 
and (d) core radius R c . 

Table 2. Independent best estimates of model parameters for terrestrial solar system planets. The preliminary reference earth model (PREM) refers 
tojDziewonski & Anderson] (| 198 l|l. 


planet 

SiQ [wt%] 

Fe/SiQ 

M g/SiQ 

Rc [km] 


reference 



Earth 

22 

0.17 

0.83 

3400 


PREM 



Moon 

23 

0.4-0.47 

0.8 

310-320 

Kuskov & Kronrod 

(2001 

) 

Mars 

21 

0.64 

0.97 

1680 

Khan & Connolly 

C 

C 

oc 


Venus 

- 

- 

- 

3100 


Aitta (2012) 



Mercury 

- 

- 

- 

2020 

Padovan et al. 

(2014) 

Hauck et al. 

(2013 



This, however, would imply a higher mantle Si-content (SIq) 

in order to simultaneously match the Fe/SiQ constraint, but 
would move SiQ outside the prior range. 

In summary, the addition of Fe/SiQ as constraint (case B) is 
clearly the key parameter that allows us to: ( 1 ) constrain man¬ 
tle composition, ( 2 ) obtain an improved estimate of core size, 
and (3) significantly reduce model variability. This ultimately 
arises because the compositional constraint (Fe/SiQ) results in 
a strong correlation of the model parameters SiQ, Fe/SiQ and 

R c in the inversion as clearly demonstrated in Figure[5] The data 
are not able to constrain Mg/SiQ, since Mg/SiQ appears not 

to be significantly correlated to the other model parameters. We 
find that independent estimates of SiQ, Fe/SiQ and R c are rel¬ 
atively well predicted by our method for Earth, Moon, Mars, and 
Venus. This is not the case for Mercury, where no model is found 
that fits data given our model assumptions. 


4. Results 


4.1. Application to synthetic exoplanets 

In this section, we apply our method to a set of synthetic 
planets where mass and radius uncertainties are considered to 
be artificially small in order to demonstrate its performance. 
We consider a range of synthetic exoplanets between 0.1 Me < 
M < 10 Me and 0.5 Re < R < 1.7 Re for both case A (in¬ 
version of mass and radius) and B (inversion of mass, radius, 
and bulk abundance constraints). Prior model ranges are listed 
in Table [T] We employ the same standard deviations used in the 
case of the terrestrial solar system planets, i.e., <jm — 0 . 01 M, 
<t r = 0, and Fe/Si 0 = Fe/Si* = 1.69 ± 0.18 (|Lodders|2003|, 
as well as fixing the minor components (Na 20 , Ca(J, and AI 2 O 3 ) 
to chondritic values (jLodders|2003|. 


Figure [TO] shows the range of synthetic rocky planets that fit 
data for case A (left panels) and B (right panels). Note that for 
every synthetic planet we solved an inverse problem, i.e., for ev¬ 
ery point in the A7-/t ) -plane in Figure]]!)] a full inverse problem is 
solved yielding separate pdfs for all parameters. When inverting 
only M and R (case A), all planets are seen to fall within the ex¬ 
treme limits of a pure iron and a pure silicate sphere (Figure [TO] 
as expected. By adding the compositional constraint (Fe/SiQ) 
inferred from observations of the host star (here the Sun), the ex¬ 
treme case of a pure iron sphere, for example, is clearly unable 
to satisfy this constraint. As a consequence, for case B synthetic 
planets plot within a smaller range than for case A. Mercury, 
for example, lies outside the permissible region for reasons dis¬ 
cussed above. 


With our assumption that planets are composed of a pure 
iron core surrounded by a rocky mantle of arbitrary size, it is 
not surprising that internal structures could only be derived for 
planets with structures lying between a pure iron and a pure 
rocky sphere. By adding a compositional constraint (Fe/SiQ) 
the range of possible internal structures is sharply reduced, lead¬ 
ing to exclusion of Mercury-type planets. 

Between the limits of a pure iron and pure silicate sphere, 
bulk densities decrease towards the silicate-sphere limit. Thus, 
we also expect core size to decrease towards the silicate-sphere 
limit. To verify that this trend is indeed obtained here, we plot 
inverted core mass fraction M c /M in Figure 11 for all consid¬ 
ered synthetic exoplanet cases shown in Figure TO In the limit of 
a pure iron (silicate) sphere M c /M —> 1 (0). For case B, models 
with M c /M > 0.6, i.e., large iron cores, are evidently excluded 
as it becomes difficult to increase mantle Si-content in order to 
satisfy the compositional constraint (Fe/SiQ). 

In order to quantify the information content of data, i.e. 
the degree to which data constrain internal structure, we com¬ 
pute the Shannon entropy measure (f7 p0 sterior/-Hprior) (Taran- 


tola|2005), which allows comparison of prior and posterior pdfs. 
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Fig. 4. Sampled ID marginal posterior distributions for Earth, (a-d) case A, (e-h) case B: prior (red), posterior (blue) and independent estimates 
(green) of model parameters: (a,e) mantle Si-content Sigj, (b,f) mantle Fe/SiQ, (c,g) mantle Mg/Si^, and (d,h) core radius R c . Independent 
estimates are listed in table[2] Note that prior ranges are different for the two cases, but axis ranges are kept the same to ease comparison. 


If .^posterior/#prior -A 1, data are only weakly constraining 
models, whereas if f? p0 sterior/-ffprior -A 0 , data pose strong 
constraints on the models (see Appendix [D] for details). Fig¬ 
ure [l 2 l shows -ffposterior/f?prior for all inverted model parame¬ 
ters. f or case A terrestrial-type planets, we find that core radius 
R c is generally well-constrained, whereas mantle composition is 
weakly constrained. For SiQ, Fe/SiQ, and Mg/Sig, the rel¬ 
ative Shannon entropy measures are larger than 0.8 throughout 
the tested range of synthetic planets. There are few planets that 
have bulk densities that are very close to the pure silicate sphere 
for which Fe/Si q appears to be better constrained. Generally 

for case A, the inversion scheme is not sensitive to variations in 
mantle composition but mostly sensitive to variations in core ra¬ 
dius R c . Stated differently, any variations in mantle composition 
can be compensated by tiny variations in core radius, while vari¬ 
ations in core radius must be compensated by large changes in 
mantle composition in order that data (M and R) be matched. 

From Figure [l2| we observe the following trend for case A: 
R c appears to be better constrained for smaller and denser plan¬ 
ets. The denser a rocky planet is, the more it is dominated by 
its core and the less mantle material is available whose compo¬ 
sitional changes could compensate for variations in core radius. 
Similarly, the smaller a planet is, the less absolute mantle mass 
is available to compensate for a higher variability in core radius. 
Therefore, the precision with which core radius can be deter¬ 
mined rests ultimately on the available mantle mass. 

For case B, the addition of Fe/Sio significantly reduces 
model variability, i.e. correlation between model parameters is 
increased. R c , SiQ and Fe/SiQ are better constrained com¬ 
pared with case A, while Mg/Sig) is only weakly constrained 
as for case A. Note however, that the prior bounds Mg/SiQ are 
significantly narrower for case B than for case A, because of the 

Article number, page 8 of 


stellar abundance constraint. The relative Shannon entropy mea¬ 
sures for core radius R c in case B are similar but slightly smaller 
than for case A. 

As discussed for the terrestrial solar system planets, only the 
core size can be constrained in the case of exoplanets for which 
only mass and radius are known (case A). In case B, the com¬ 
positional constraint (Fe/SiQ) introduces a strong correlation 
between the model parameters SiQ, Fe/SiQ and R c . This en¬ 
ables us to constrain mantle composition. 


4.2. Application to confirmed exoplanets 


In this section we apply our methodology to actual exoplanet 
measurements in order to derive models of the internal structure. 
At this point in the development of our approach we are still lim¬ 
ited by the assumptions made for the overall structure (iron core 
and rocky mantle with no ice and no atmosphere) and hence this 
application serves more as an illustration of the overall capabili¬ 
ties of our method. In Figure[l3]we plot all confirmed exoplanets 
in the mass-radius range of interest on top of our synthetic ex¬ 
oplanet calculations (Figure TO]). Masses, radii, and uncertainty 
ranges of the observed exop anets are listed in Table [3] Kepler- 
680 , Kepler-131c and Kepler-406c have estimated mean densi¬ 
ties above that of pure iron and clearly present special exotic 
cases. Kepler-36b has been chosen to demonstrate the applica¬ 
tion of our method to actual observations because of small un¬ 
certainties on mass and radius. To our knowledge, no stellar ele¬ 
mental abundances of Fe/SL*. and Mg/SL*. of its host star other 
than for the Sun are available. Hence, we test case B assuming 
solar elemental abundances. Since the Kepler mission has been 
targeting Sun-like stars (Borucki et al .|2010| >, this is a permissi¬ 
ble assumption. 
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Figure [l4| shows the results for Kepler-36b. For case A, we 
see that the mantle ratios Fe/SiQ and Mg/SiQ are not con¬ 
strained at all, while SiQ and R c are moderately and well con¬ 
strained, respectively. By adding Fe/Sio to the data (case B), 
SiQ and Fe/SiQ are significantly better constrained. Thus, in¬ 
corporating observations of stellar elemental abundances has a 
big impact on our ability to constrain model parameters. The 
posterior range of R c is very large in both cases ranging from 
zero to about half the planet radius. We argue that mass and ra¬ 
dius uncertainties of 7 % and 2 %, respectively, are the reason 
for such a large spread. 

The example of Kepler-36b illustrates the necessity of having 
well-measured masses and radii in order to derive constraints on 
the interior structure. Furthermore, it exemplifies the large im¬ 
pact of elemental abundance ratios in order to constrain plane¬ 
tary interior structure. Since it will be impossible to directly mea¬ 
sure these, spectroscopically-determined values measured from 
the photosphere of the host star should serve as proxies. 


4.2.1. Information content of the data 

The limiting factor on our ability to constrain interior struc¬ 
ture is measurement precision. In this section, we study in more 
detail the link between measurement precision and resulting con¬ 
straints on internal structure. 

For this, we compute the relative Shannon entropy measure 
^posterior/-ffprior for each model parameter for a range of mea¬ 
surement uncertainty. In these calculations, inversions were per¬ 
formed using a grid-search method. The results are shown in 
Figure [15] for the example of Kepler-36b and in Figure [16] for 
an Earth-like planet. The results can be summarized as follows: 


- Overall, a decrease in data uncertainty leads to a smaller 
/^posterior//7prior and thus better constrained parameters. 
However, for large uncertainties H poster i or /H plIOI is nearly 
constant and relatively large. 

- Uncertainties on M, R, and Fe/SiQ appear equally impor¬ 
tant. 

- The inherent degeneracy of the problem is limiting our 
ability to constrain interior structure even in the case 
of very small data uncertainties. For the cases shown, 
-ffposterior/-ffprior does not fall below 0.5. 

- Our ability to constrain interior structure is significantly dif¬ 
ferent for larger and denser exoplanets such as Kepler-36b 
compared to Earth-like bodies. As shown for Kepler-36b, 
fTp OS terior/-ffprior decreases only at extremely small data un¬ 
certainties whereas for an Earth-like body, -ff pos terior/-ffprior 
decreases significantly earlier. Given their large difference 
in mass, this indicates that measurements with a given pre¬ 
cision are able to better constrain the internal structure of 
small-mass planets than large-mass planets. Large planets 
have more absolute mantle mass available whose composi¬ 
tional changes are able to compensate for a higher variabil¬ 
ity in core radius, ultimately reducing our ability to constrain 
interior models. 

- Over the entire tested range of uncertainties on Fe/SiQ (0- 
15%), significant improvements in our ability to constrain in¬ 
terior structure can be observed, whereas for mass and radius 
this is only evident in the case of small uncertainty ranges 
(< 5%). 


5. Discussion 


The results presented here have to be interpreted in the light 
of the assumptions made. Firstly, we consider terrestrial-like 
planets that consist only of silicate mantles and pure iron cores. 
Sotin et al. ( 2007| showed that rocky planets can be described 
with great accuracy by only using the four elements Si, Mg, Fe, 
O. Secondly, we assume a pure iron composition for the core. 
Volatiles (H, He, H 2 0, etc.) are currently neglected. Compared 
to Earth, where the core include lighter elements (e.g., S, Si, 
C, O), our method therefore systematically overestimates core 
density and underestimates core radius. This introduces biases 
in the estimates of mantle composition due to the strong cor¬ 
relations between mantle composition and core radius. In the 
future, lighter elements (e.g., Si) will be included in the core. 
Thirdly, our approach assumes subsolidus conditions and per¬ 
fectly known Equation-of-State parameters for all considered 
compositions. For exoplanets that are larger than Earth, extrapo¬ 
lations of parameters to pressures and temperatures beyond those 
measured in the laboratory are required. This introduces uncer¬ 
tainty in the EoS modeling that we will address in the future. 

While the mean density of a terrestrial-like planets can be de¬ 
rived from measurements of mass and radius, its internal struc¬ 
ture is only relatively poorly constrained by these two quantities 
as also shown by Rogers & Se ager| (e.g., |2010| ). We find that 
adding Fe/SiQ and Mg/SiQ as compositional constraints sig¬ 
nificantly improves this situation. These additional constraints 
significantly narrow the prior range of mantle compositions (for 
details on use of different priors see e.g., Rogers & Seager :2010| l. 
In particular, Fe/SiQ introduces a strong correlation between 
model parameters SiQ, Fe/SiQ, and R c , which significantly 

reduces the range of possible models. Hence, Fe/SiQ plays a 
key role in constraining mantle composition and internal struc¬ 
ture. 

Our ability to constrain interior structure is directly limited 
by data uncertainties. We have demonstrated that uncertainties 
on mass, radius, and elemental abundances appear to be equally 
important. While it is worth striving to increase the measure¬ 
ment precision on all input data, the knowledge of Fe/Si^ and 
Mg/SL*. is particularly important as it represents a relatively 
powerful means of reducing the degeneracy of possible interior 
models. The level of inherent degeneracy in model solutions ac¬ 
tually depends on the specific values of M and R (as discussed 
for Figure [T2| ( [Rogers & Seager|2010| : for given measurement 
accuracy of mass and radius, interior structure is better con¬ 
strained for an Earth-sized planet than for a larger body (e.g., 
Kepler-36b). 

In an earlier study, [Rogers & Seager] ( ]2010| employed 
Bayesian analysis to study exoplanet GJ 581d in an attempt 
to constrain core mass and silicate mantle mass fractions of 
MgSiOs and FeSiOa. Based on mass and radius alone, [Rogers] 
& Seager ( |2010] ) also concluded that an exact interior solu¬ 
tion cannot be inferred since the problem is highly underde¬ 
termined. However, Rogers & Seager] ( ]2010| l argue against the 
use of bulk abundance constraints inferred from Fe/SL*. and 
Mg/SL*. in order to remain as independent as possible of planet 
formation models. It should be noted that if Rogers & Seager 
( |20T0| were to use bulk abundance constraints core size and 
mantle composition would be fixed since they use a parameter¬ 
ized phase-diagram approach and restrict the mantle mineralogy 
to Mg ] _ Y Fe A SiO ;i . Consequently, the problem would no longer 
be degenerate. Clearly, our approach is more general because it 
allows for arbitrary elemental abundance ratios while using stel- 
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Table 3. Masses and radii of observed exoplanets 


planet 

radius R/Re 

mass M/Me 

reference 


mean bulk density g/cm 3 

CoRoT-7b 

1.55 ±0.10 

7.31 ±1.21 

Moutou et al. 

2013; 


10.82 

Kepler-10b 

1.416 ±0.033 

4-5 QtlH 

Batalha et al. ( 

2013) 


8.89 

Kepler-36b 

1.486 ±0.035 

4.4518:2? 

Carter et al. (. 

> 012 ) 


7.48 

Kepler-5 7 c 

1.55 ±0.04 

5.4 ±3.7 

Steffen et al. 

2013) 


8.00 

Kepler- 68 c 

0.953±g;8S 

4-81^ 

Gilliland et al. 

(p013 

) 

30.58 

Kepler-7 8 b 

1 1 '7Q+0-159 

± * i ' °_0.089 

i-86l8:2s 

Pepe et al. (2 

013) 


6.35 

Kepler-93b 

1.50 ±0.03 

2.59 ± 2.00 

Weiss & Marcy 

(2014 

b 

4.23 

Kepler-97b 

1.48 ±0.13 

3.51 ±1.90 

Weiss & Marcy 

(2014 

b 

5.97 

Kepler-99b 

1.48 ±0.08 

6.15 ±1.30 

Weiss & Marcy 

(2014 

b 

10.46 

Kepler-100b 

1.32 ±0.04 

7.34 ±3.20 

Weiss & Marcy 

(2014 

b 

17.60 

Kepler-102b 

1.18 ±0.04 

3.8 ±1.8 

Weiss & Marcy 

(2014 

b 

12.75 

Kepler-131c 

0.84 ±0.07 

8.25 ±5.90 

Weiss & Marcy 

(2014 

b 

76.74 

Kepler-406b 

1.43 ±0.03 

6.35 ±1.40 

Weiss & Marcy 

(2014 

b 

11.97 

Kepler-406c 

0.85 ±0.03 

2.71 ±0.80 

Weiss & Marcy 

(2014 

b 

24.33 


lar abundance constraints as a means of further reducing model 
variability. 

With regard to planet formation, these studies suggest that 
some interior compositions are more likely than others. Recent 
planet formation studies accounting for equilibrium condensa¬ 
tion of material (e.g., Thiabaud et al.|2014 1 have shown that most 
planets have a bulk refractory (Fe, Si, and Mg) composition that 
appears to be indistinguishable from that of the host star because 
species of Fe, Si and Mg condense at similar temperatures (~ 
1000 K). Hence, most planets will have similar refractory ele¬ 
mental ratios as those in the original disc which is assumed to 
be well-represented by stellar ratios. Also, O, Fe, Si, and Mg 
are predicted to be among the most frequent solid species in cir- 
cumstellar disks (e.g., Elser et al. |2012| |Gilli et aLp 006 ) as a 
consequence of which models with large iron cores and small 
silicate mantles are generally unlikely to form. Based on these 
observations Valencia et al. ( [2007b i also proposed the general 
use of a minimum bulk value for Fe/Si. Here, we have followed 
Grasset et al.; ( [2009) , who considered bulk Mg/Si and Fe/Si to 
be dictated by their abundance in the host star’s photosphere and 
adopt these as proxies for planet bulk composition. 

For the terrestrial solar system planets, we showed that us¬ 
ing solar photospheric abundances (Fe/Si* and Mg/Si*) yields 
tighter constraints on the internal structure, except for Mercury. 
Mercury’s interior is unique among the terrestrial planets and is 
usually explained by the removal of a large fraction of its sil¬ 
icate mantle by either a giant impact or by evaporation due to 
the proximity of the planet to the sun ( |Benz et aI7||1988) >. Con¬ 
sequently, our scheme is not able to find interior models that fit 
data for Mercury. In spite of this, we argue for the use of stellar 
abundance constraints, because of the ability of our scheme to 
predict for which planets alternative accretion scenarios or post¬ 
formation processes must be envisaged. We also note that mete¬ 
orites, which are assumed to be the building blocks of the terres¬ 
trial planets, have abundances of refractory elements within 10 % 
of solar values (]Lodders|2003|). 


6. Conclusions and Outlook 

We have presented a Bayesian inversion method in order to 
constrain the interior structure of rocky exoplanets from mea¬ 
surements of mass and radius. Since stellar refractory elemental 
abundances are likely to be generally a good proxy for plane¬ 
tary bulk compositions, they represent powerful additional con- 

Article number, page 10 of 


straints allowing us to significantly improve determination of in¬ 
ternal structure. Our proposed scheme is more general than pre¬ 
vious works on mass-radius relationship of rocky exoplanets, in 
that we obtain confidence regions of interior structures of general 
mantel composition and core structure. Furthermore, we investi¬ 
gated how interior structure models depend on various parame¬ 
ters in the form of prior information, data, and data uncertainties. 

We have applied the method to the terrestrial planets of our 
solar system and compared the results to independent estimates 
and, except for Mercury, find generally good agreement. Follow¬ 
ing this, we have applied our method to synthetic exoplanets in 
a large mass and radius range as well as to exoplanet Kepler- 
36b. We have investigated two different cases: case A where we 
invert mass and radius, and case B where we additionally con¬ 
sider stellar elemental abundances obtained from photospheric 
observations of the host star. We summarize our findings in the 
following: 

- Photospheric elemental abundances of the host star (Fe/Si* 
and Mg/Si*) are key parameters in reducing model de¬ 
generacy through the introduction of correlations between 
mantle composition and core size. 

- How well mantle composition and core radius can be gen¬ 
erally constrained depends on data and data uncertainties. 
There is an inherent degeneracy that limits our ability to 
constrain the interior structure even in the case of small data 
uncertainties. Independently of these, it seems that model 
variability depends on mass and density of a planet. Thus, it 
is our contention that a case-by-case probabilistic inversion 
that provides model parameter uncertainties is indispensable 
in order to rigorously characterize interior structure. 

- Measurement uncertainties on mass, radius, and stellar 
abundance constraints appear to be equally important. 

- For improved characterization of interior structure better es¬ 
timates of Fe/Si* and Mg/Si* are required. We have shown 
for Kepler-36b and an Earth-sized planet, that significant im¬ 
provements on model parameter estimation can be achieved 
by reducing the uncertainty of Fe/Si* and Mg/Si*. 

Space missions that aim at characterizing exoplanets by 
means of the transit method seek the most precise measurement 


19 































































Dorn et al.: Can we constrain interior structure of rocky exoplanets from mass and radius? 


of R and hope for follow-up missions to provide precise mea¬ 
surements for M. However, since data precision also depends 
on characteristics of target star and observation time (integra¬ 
tion time, number of transits observed), a careful weighing of 
costs and benefits is crucial for the success of a mission. We 
have demonstrated how the relative Shannon entropy can be used 
in future missions (e.g., CHEOPS, TESS, PLATO) as a means 
of optimizing the scientific return. In particular, we are able 
to quantitatively assess the improvement in interior structure- 
determination of a specific exoplanet target due to an increase in 
the measurement precision of mass and radius. To illustrate this 
we consider the case (discussed in Figure [l5| l, where the pre¬ 
dominantly flat curves imply that smaller uncertainties on mass 
and radius do not provide significantly better insights on the in¬ 
terior structure. For an Earth-sized exoplanet, however, smaller 
uncertainties on mass and radius can lead to an improved under¬ 
standing of interior structure (Figure [l 6 |. 


An increasing number of Super-Earths are being observed. 
Super-Earths lie in the intermediate mass-range between terres¬ 
trial planets and the gas/ice giants in the solar system with dif¬ 
ferent scenarios for their interiors. Their interior structures (e.g., 
purely rocky composition or predominantly water/carbon com¬ 
pounds) and formation histories are, however, a matter of debate, 
as is the question of whether Super-Earths can harbour life (e.g. 
Valencia et al. 2007a Howe et al. 2014|>. A potential way of 


answering these questions is to invert for the physico-chemical 
structure in order to narrow down the possible types of interior 
models. 


In future studies, we will extend the dimensionality of our 
problem to include more elements, e.g., hydrogen and water, 
so as to model also exoplanets with atmospheres and/or oceans. 
This will eventually allow us to study how interior types are dis¬ 
tributed among stars and what can be learned from such distri¬ 
butions about planet formation. However, for a given data set, 
increased model complexity is generally accompanied with a 
larger degeneracy in model solutions. Additional data would 
help to restrict model variability. For example, dynamic obser¬ 
vations can constrain a planet’s tidal dissipation and thereby 
provide constraints on the distribution of rigidity and density 
(through the Love number, e.g.. Delisle et al.|2014l. 
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Appendix A: Prior of mantle Si-content 

The prior bounds on Sig depend on the prior bounds of 
Fe/SiQ and Mg/Si^, since the mass fractions of all oxides 

in the NCFMAS chemical system must sum to one. The mass 
fraction of Si can be expressed as 


Si 


0 we find the prior of Si^ being 


P( Si @) I 


<9Fe/Si@ 
9Si@ 1 


(1 ~ XNa 2 0,Ca0,Al 2 0 3 ) 



A'FcO 

MFe 


This implies a non-uniform prior for Si^. Outside the given 
interval in Table [Tj the prior distribution of Si^ is essentially 
zero. 1-1 


Appendix B: Temperature profile 

We assume a fixed mantle temperature profile that is based 
on the Earth model. A variable temperature profile for the mantle 
only introduces negligible variations in the density profile (< 1 
%). This is demonstrated in Figure |B. 1 1 Surface temperature is 
fixed at 1000 K and although temperature varies between 2100 
K to 3500 K at the core mantle boundary, resultant changes in 
density are relatively small (Figure |BTT|>. 


Appendix C: Silicate model chemical system 

To illustrate the influence of the chosen model chemical sys¬ 
tem, we solve the inverse problem for the case of the Earth as in 
section 3.2, but using the simpler FMS system. Sampled model 
parameter pdfs are shown in Figure |CT| and relative to NCFMAS 
(Figure[4]), we observe that with the FMS system mantle compo¬ 
sition and core radius are less well-predicted, although poste¬ 
rior model parameter variability is, as expected, smaller (fewer 
chemical components). Also, future use of stellar determinations 
of Ca, Na, and Al abundances support the use of the NCFMAS 
system. 


Appendix D: Shannon entropy measure 

The Shannon entropy measure of a discrete variable X of 
possible states {x\, X 2 , Xn}, with N being the number of 
bins, can be written as: 


N 

H = ~J2p(Xr) * log 2P(Xi) 

z — 1 

where the Xi states represent the distribution histogram bins in 
which the probability of occurrence is given by 

^ # of realizations in the bin 

1 total # of realizations 

The Shannon entropy H has a maximum fT max = log 2 (-/V) 
when all N states are equally likely (e.g., a uniform distribution) 
and a minimum H min = 0 when all realizations fall into a 
single histogram bin. Clearly, H is bin-size dependent. Here, the 
bin size used is defined such that the prior range of each model 
parameter is divided in « 20 different bins. 


Si ® 


( 


MSiQ 2 

MSi 


(1 - XNa 2 Q,Ca0,Al 2 0 3 ) _ 

Fe / Si ®^+ M g/ Si ®^) 


where // are respective molar weights and the mass fraction of 
minor oxides are combined in XNa 2 0,Ca0.Ai 2 0 3 - Transforming 
from a uniform prior in Fe/Si^ and Mg/Si^ to the variable 
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Fig. 5. Sampled 2D marginal posterior distributions for the Earth (case Agplots a-f and case B: plots g-1) showing correlation between parameters 
Sig), Fe/Sijg), Mg/Sig), and core radius R c (corresponding to Figurepj. Independent estimates (TablepJ are marked in white. 
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15 20 25 0 1 2 0.7 0.8 0.9 1 1.1 0 500 1000 1500 


mantle Si [wt%] mantle Fe/Si mantle Mg/Si core radius R c [km] 

Fig. 6. Sampled ID marginal posterior distributions for the Moon, case B: prior (red), posterior (blue) and independent estimates (green) of model 
parameters: (a) mantle Si-content Si^, (b) mantle Fe/Sigj, (c) mantle Mg/Sijg,, and (d) core radius R c . Independent estimates are listed in 
table [2] 
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Fig. 7. Sampled ID marginal posterior distributions for Mars, case B: prior (red), posterior (blue) and independent estimates (green) of model 
parameters: (a) mantle Si-content Sigj, (b) mantle Fe/Sigj, (c) mantle Mg/Sigj, and (d) core radius R c . Independent estimates are listed in 
table [2] 



Fig. 8 . Sampled ID marginal posterior distributions for Venus, case B: prior (red), posterior (blue) and independent estimates (green) of model 
parameters: (a) mantle Si-content Sig), (b) mantle Fe/Sigj, (c) mantle Mg/Sigj, and (d) core radius R c . Independent estimates are listed in 
table [2] 
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Fig. 9. Sampled ID marginal posterior distributions for Mercury, case B: prior (red), posterior (blue) and independent estimates (green) of model 
parameters: (a) mantle Si-content Sigj, (b) mantle Fe/Sigj, (c) mantle Mg/Sig), and (d) core radius R c . No interior model fits the data, i.e., 
likelihoods are nearly zero. Independent estimates are listed in table [2] 
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Fig. 10. Range of synthetic rocky planets that fit data (black dots) of (a) M and R (case A), and (b) M, R, and stellar abundance constraints 
(case B). Mass-radius curves for pure iron (red) and pure silicate spheres (blue, zero amount of FeO) were calculated independently. In case A, all 
synthetic planets that fit data (black dots) fall within the limits of a pure silicate (blue line) and pure iron sphere (red line). In case B, this region is 
reduced to a smaller area indicating that chemical bulk ratio (Fe/Sio) introduces strong additional constraints on the internal structure. Synthetic 
planets for which no internal structure could be derived within model assumptions (pure iron core surrounded by a rocky mantle) are indicated by 
gray open circles. Note that for every circle in the plots an inverse problem has been solved. See main text for further details. 
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Fig. 11. Mean posterior core mass fraction M c /M plotted for the range of synthetic rocky planets that fit (a) M and R (case A), and (b) M, R, 
and stellar abundance constraints (case B). (See also Figure [T()|) 
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Fig. 12. Relative Shannon entropy measure //posterior/-//prior of model parameters plotted for the range of synthetic rocky planets that fit a given 
Aftfiild A^Sgfland M, R, and Fe/Sio (case B, right panels). Model parameters are well constrained for //posterior///prior —> 0 

(blue colors) and not constrained for //posterior///prior —> 1 (red colors). For case B, variations of Z/ P osterior/Z/prior along the grid appear not to 
be smooth which is due to the poor grid resolution. See also Figure [Hi] 
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Fig. 13. Confirmed exoplanets plotted against Figure[lO] (a) Case A and (b) case B. 





Fig. 14. Sampled ID marginal posterior distributions for Kepler-36b, case A (a-d) and B (e-h): prior (red) and posterior (blue) of model parameters: 
(a,e) mantle Si-content Sigj, (b,f) mantle Fe/Sigj, (c,g) mantle Mg/Sig), and (d,h) core radius R c . Data uncertainties are <j Fe /g io = 10%, 

<jr = 2%, and vm = 7%. 
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Fig. 15. Summary of information content for Kepler 36-b (M = 4.45Mb, R = 1.486-Rb): Relative Shannon measure as a function of (a) a m 
while fixing or at 0.5% (b) cr Fe /Si 0 while fixing <tm at 0.5% (c) or while fixing <TFe/Si 0 at 0.5%. 



Fig. 16. Summary of information content for an Earth-sized body (M = 1 Me, R = 1 Re): Relative Shannon measure as a function of (a) um 
while fixing ctr at 0.5% (b) a-Fe/Si 0 while fixing Om at 0.5% (c) ctr while fixing <TFe/Sio at 0-5%. 
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depth [km] 

Fig. B.l. Sampled profiles of density profiles for Earth. A variable mantle temperature (gray profiles) only introduces little additional variation in 
density compared to profiles with fixed mantle temperature (green colors). 



Fig. C.l. Sampled ID marginal posterior distributions for Earth using the FMS system: prior (red), posterior (blue) and independent estimates 
(green) of model parameters: (a) mantle Si-content Sigj, (b) mantle Fe/Sig), (c) mantle Mg/Sigj, and (d) core radius R c . Independent estimates 
are listed in table [2] 
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